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ABSTRACT 

Using the reduced formulation of large- A Quantum Field Theories we study 
strings in space-time dimensions higher than one. We present results on possi- 
ble string susceptibilities, macroscopic loop operators, l/A^ 2 -corrections and other 
general properties of the model. 
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1. INTRODUCTION AND CONCLUSIONS 



One of the most interesting open problems in String Theory and Quantum 
Field Theory is the study of the systems with string excitations in dimensions 
below the critical dimension. In the context of lattice gauge theories the problem 
can be formulated in terms of the Schwinger-Dyson equations satisfied by the 
Wilson-loop observables. In the context of large- N field theory, it was noticed in 
the early eighties that due to the factorization properties of most gauge invariant 
observables, the loop-equations of the theory were equivalent to those of a theory 
on a single hypercube with periodic boundary conditions. These models are known 
as the large- N reduced models [1,2,3]. 

In this paper we use the techniques of the reduced models to study non-critical 
strings in dimensions higher than one. We will find an action in terms of a single 
matrix which exactly reproduces the partition function of a string moving in 2, 3, 4 
space-time dimensions. The model does not appear yet to be exactly solvable, 
and we can at the moment make only general preliminary remarks concerning its 
critical properties and its possible continuum limit. 

The first lattice formulations of string theories in arbitrary dimensions ap- 
peared in the mid-eigthies [4,5,6,7], and there was a good deal of activity in 
solving some two-dimensional models on random triangulated surfaces [10]. A 
breakthrough took place with the work in ref. [8] which provided a description in 
the continuum of the coupling of minimal conformal models to two-dimensional 
gravity in the light-cone gauge. Their subsequent formulation in the conformal 
gauge, and the generalization of the results to surfaces of arbitrary topology ap- 
peared in [9]. With the discovery of the double scaling limit [11] we begun to un- 
derstand some of the perturbative and non-pertubative properties of non-critical 
strings in dimensions below or equal to one. Little progress has been made however 
in going beyond d — 1 apart from some evidence that it is possible to construct 
Liouville Quantum Field Theory in some special dimensions d = 7, 13, 19 . Some 
interesting recent results appear in [13]. 
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Some interesting problems involving the presence of external field have also 
been studied [14,15,16,17]. In particular the work in [17] attempted to generalize 
the standard one-matrix model to include the presence of local curvatures in the 
triangulated surfaces described by the large- N limit. We will make contact with 
this work later in sections three and four. 

The outline of this paper is as follows: In section two we formulate the problem 
of strings in arbitrary dimensions in general, we present the general methodology 
of reduced actions, and the explicit form of the one-matrix model action which 
reproduces exactly the properties of planar string in .D- dimensions. Some details 
are given concerning the way to carry out perturbative computations in the context 
of reduced models. 

In section three we continue the analysis of our model and show how in the 
simplest possible approximation one obtains a second order phase transition be- 
tween a pure gravity phase with string susceptibility 7 s t = —1/2 and a branched 
polymer phase with 7 s t = +1/2. This is in agreement with the analysis carried out 
in [17] for the one- matrix model with world-sheet curvature. 

In section four we analyze some general effective actions which could be taken 
as approximations to our theory. We solve the planar limit of these theories exactly 
by looking for one-cut solutions. We find generically a second order phase transition 
between a pure gravity-like phase with negative 7 s t (with the value 7 s t = — l/(m + 
1), as in the Kazakov multicritical points [7], and a phase that resembles branched 
polymers with positive string susceptibility 7 s t = n/(n+m+l), with n, m arbitrary 
positive integers. In the approximation we envisage, the D-dimensional embedding 
of the string is effectively replaced by some kind of polymeric matter whose role is 
to break the world-sheet into surfaces oscullating at several points. We have the 
following phases: a phase where gravity becomes critical and matter remains away 
from criticality, in which case we obtain the pure gravity regime; a phase where 
matter becomes critical while gravity remains non-critical, leading to a theory 
of branched polymers; and finally the most interesting phase where both become 
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critical simultaneously. This phase is described by the positive string susceptibility 

7st = n/(n + m + 1). 

To understand better the geometrical properties of this phase, we compute 
correlators of macroscopic loop operators on the surface in section five. We show 
that it is possible to define macroscopic loops in the continuum limit only for a 
subclass of critical points with n — 1. We find it quite remarkable that only for 
the n = 1 critical points is there enough room on the surface to open arbitrary 
numbers of loops of arbitrary macroscopic length. For higher values of n, the 
polymerization of the surface is presumably so complete that no room is left to 
construct such operators. This leads us to believe that only in this phase we may 
be able to recover some of the basic properties of non-critical strings. 

The continuum limit for n — 1 contains an extra state with respect to the 
standard pure gravity case, representing the breaking of the surface into two pieces 
touching at one point. This extra state is probably the shadow of the effect of the 
tachyon that one expect to find as soon as string theories are formulated in dimen- 
sions higher than one. Whether this feature remains when more realistic effective 
actions are considered remains to be seen, but we find it rather encouraging. We 
also explore some of the scaling operators appearing in truncations of the reduced 
model. 

To make sure that the scaling limit we take indeed defines a string theory, in 
section six we solve the Schwinger-Dyson equations for the approximate effective 
action on genus one. We find that the David-Distler-Kawai dependence of the 
string susceptibility on the genus of the surface is recovered to this order. This is 
encouraging evidence that the continuum limit that we define does describe some 
kind of non-critical strings. 

Although we find the results just summarized encouraging, there are several 
obvious points in our approximations which should be clarified. First of all, in the 
effective actions studied in section four there is little evidence of the dependence of 
the string susceptibility on the string embedding dimensions. In particular, in the 
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reduced model one of the couplings represents the two-dimensional cosmological 
constant. There is a similar coupling in the effective action for the reduced model 
which plays a similar role. What needs to be done is to show that the relation be- 
tween these two couplings is analytic, at least for some range of parameters. The 
conclusions we draw in sections four and five depend implicitly on this assumption, 
and this is an important caveat to keep in mind when considering the work pre- 
sented. The most natural thing to do seems to be to try to settle these questions 
by a direct renormalization group analysis of the reduced model along similar lines 
to those suggested in [18], or to construct an approximate Migdal-Kadanoff-type 
transformation. Work in this direction is in progress. We hope to convince the 
reader that the study of non-critical strings in dimensions higher than one, using 
reduced large- N model is worth pursuing. 

While this work was being completed we received several preprints having some 
overlap with our work [19]. 

2. LARGE-TV REDUCED THEORIES AND MATRIX MODELS 

In the study of strings propagating in flat D-dimensional space-time, we repre- 
sent the sum over two-dimensional metrics in Polyakov's approach to string theory 
[20] by a sum over tringulated surfaces. The quantity we would like to evaluate 
[4,6,7] is 




Where g is the genus of the triangulation. For a given triangulation T, Tb,Ti,T2 
are respectively the sets of vertices, edges and faces of T, n(T) is the order of the 
symmetry group of T, |T| is the total area of T counting that every triangle has 
unit area. \i is the bare two-dimensional cosmological constant, X{ describes the 
embedding of the triangulation into D-dimensional flat space-time and G(X — Y) is 
the propagator factor for each link. We have included also the local volume factor 
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of at site % to represent the effect of local curvature. Usually we take a = D/2, 
but it is more convenient to leave this exponent arbitrary to include the effect of 
local curvature terms on the world-sheet. For a triangulation, if qi is the local 
coordination number at site % (the number of triangles sharing this vertex), 



Ui = 3' 



and the volume factor in the measure can be shown to be related in the continuum 
to an expansion of the form: 



^2 m °i = c o + ci J JgR + c 2 J ^/gR 2 



thus representing some of the effects of world-sheet curvature. The general features 
of the phase diagram in the (D,a) plane were studied for example in [21]. In the 
standard Polyakov formulation, the propagator is Gaussian, 

G(X) = e-* 2 /2 

By standard large-N analysis the sum (2.1) can be transformed [6,7] into the large- 
N expansion of a matrix field theory in D-dimensions 

Z= lim log / Z^expf-iV / d D XTr(-(j)G~ 1 (j) + -# 3 ) J (2.2) 
N->oo J \ J 2 3 / 

It is difficult to proceed very far with Gaussian propagators. If we work in dimen- 
sions D < 6 it should not matter whether we replace the Gaussian propagator by 
the Feynman propagator. There is numerical evidence that this change does not 
affect the critical properties of the theory below six dimensions [22]. Notice that 
in this construction the exponent a in (2.1) is set to zero. If we want to include 
the effect of local world-sheet curvature, we can follow [17] and change the kinetic 
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term in (2.2). The two-dimensional cosmological constant is here represented by 
g, g = e~ M . To summarize, we want to study the critical properties of the action: 

Z = J Dtf>(X) exp (-N J d D XTr( 1 -d fI <j ) d^ + |0 3 )) , (2.3) 



or 



' ;! ! ) , (2-4) 



Z = J D(f){X) exp(^-N j d D XTr( y ^Ad )I (pAd fl (p + | 

where A is a constant N x N matrix and is an iV x A" matrix field. The reason 
why the A-matrix in the kinetic term simulates the effect of local curvature can be 
seen by writing the propagator in (2.4) explicitly 

{MXWuOO) = AjiA^G(X - Y). 

If we ignore for the time being the propagator factor G, for every closed index loop 
in a generic 3 graph made of q propagators, we obtain a contribution of TrA~ q . 
Since the 3 graphs are dual to triangulation, this means that we are associating a 
curvature factor of TrA~ q to the vertex dual to the face considered. In this way we 
can simulate the presence of the a a term in the measure in (2.2). We are interested 
in particular in the computation of the string susceptibility exponent 



dp 



x~^r -(g-gc)-^, (2.5) 



where F is the free energy of the system, and g c is the critical value of g indicating 
the location of the critical point. Apart from the string susceptibility we would 
also like to compute the spectrum of scaling operators and the properties of the 
quantum geometry implied by (2.3) (2.4). 

To simplify the arguments we concentrate on the study of planar configurations 
(spherical topologies). From work dating back to the late seventies [23], and due 
to the factorization properties of the leading large-A^ approximation, the planar 
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limit is dominated by a single constant configuration; Witten's master field (a 
master orbit in the case of gauge theories). This idea is explicitly realized in 0- 
dimensional matrix models [24]. In the lattice gauge theories one can think of 
the reduced Eguchi-Kawai (EK)[1] or Twisted Eguchi-Kawai (TEK)[2] models as 
explicit descriptions of the master field. We now briefly describe the main ideas 
behind the reduced models. Further details can be found in the previous references 
and in the review article [3]. These models were originally formulated to describe 
the large- N properties of lattice gauge theories. 

Due to the factorization property in the large- N limit, the Schwinger- Dyson 
equations for Wilson- loops form a closed system of equations in the planar limit. 
For gauge theories we get an infinite set of polynomial equations for the Wilson- 
loops [25]. In terms of the standard link variables U^{x), the EK proposal consists 
of making them independent of x, U^(x) i— > U^, and the action of the reduced 
model becomes 

s EK = ^S(U ll (x)~U ll ). (2.6) 

Some of the properties of this action are: 

1) . We obtain a theory on a single hypercube with periodic boundary condi- 
tions. This is an important simplification of the problem. 

2) . In the case of gauge theories, the gauge symmetry becomes a global sym- 
metry, Ufj, i— > QU^Q -1 . 

3) . There is an extra U(1) D symmetry in the case of U (N) lattice gauge theory 
Up i-> e^Up. In the case of SU(N) the symmetry becomes Zfy. 

4) . It is possible to show that the planar Schwinger- Dyson equations following 
from the reduced action coincide with those obtained from the original Wilson 
theory as long as open loop expectation values vanish. This can be shown to 
be true at strong coupling, but it does not hold at weak coupling [26]. Without 
this problem we would have a beautiful implementation of the Master Field idea, 
because the loop equations for the 1-site EK model are identical with the standard 
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ones. The problem is connected with the fact that at weak coupling the extra 
U(1) D symmetry is broken. Open loops have non-trivial charge with respect to 
this symmetry. Only if the symmetry remains unbroken are we guaranteed to 
maintain the loop equations without extraneous terms. 

The resolution of this problem motivated the formulation of the TEK model 
[2]. It is inspired on 'tHooft's use of twisted boundary conditions in gauge theories 
[27]. For any matrix theory the TEK prescription is exceedingly simple. We reduce 
according to 

4>{x) h-> D{x)<j>D{x)- 1 , (2.7) 

where D(x) is a projective representation of the D-dimensional lattice translation 
group. We have to choose a set of D iV-dimensional matrices such that: 

D(x)=Y[l%' 1 , x=( Xl ,X2,...,x D ), Xi eZ, (2.8) 

since D(x) has an adjoint action on the fields 0, the T^'s are required to commute 
only up to an element of the center of SU (N) , 

r^IV — ZvjjTjsT ^, 

= e 2nin ^/ N , (2.9) 
and the integers n^ v are defined mod N . The reduced action prescription is now 

Stek = ^[ S ^( x ) ^ Dix^Dix)- 1 ), (2.10) 

and similarly for expectation values. 

To avoid open loop expectation values and a mismatch between the original and 
reduced Schwinger-Dyson equations, the matrices T M must verify some conditions. 
In particular they should generate an irreducible representation of the group of 
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lattice translations. For a lattice with LP sites this requires N = L D I 2 . Hence 
N = L for D = 2 and N = I? for D = 4. Hence, if we want to simulate a 
two-dimensional lattice with 64 x 64 sites it suffices to consider the group SU(64). 
The choice of twist matrix depends on the dimensionality. In two-dimensions, 
the simplest choice is given by 

n^y = e^y. (2.11) 

The explicit form of the four and higher dimensional twists can be found in the 
quoted literature. 

We now set aside lattice gauge theory and return to our problem. The twisted 
reduced version of (2.3) becomes: 

S = ^Trp^r- 1 - 4>? + TrV{4>). (2.12) 
n 

The equivalence of the planar approximation to (2.3) with (2.12) follows from ear- 
lier work on reduced models [2,28]. For simplicity we consider the two-dimensional 
case. First, since for = e^qy and q v defined modulo N, 

A(q) = T^ 2 

are N 2 linearly independent matrices, we can expand the matrix in the T-basis 
according to 

= ^(f) q A{q). 
Some useful properties of the A-matrices are: 

A((/)t = A(-q)e™(k\k), (kW) = J2 n ^ k i = k 2 k i 

A(q\ ...,q n )= A{q l + • • • + q n ) exp I ^ J>W> 

\ i<3 
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and 

ryi^rt = e 2 ^ L A( q ), 



explicitly showing how the T's implement the lattice translation group. Further- 
more, 

Tr(A(qU( q '))=N5 L (q-q'), 



TriA(q 1 )...A(q n )) = N6(J2^^ £< W j > 

where the 5-function is defined mod L. Using these properties, the kinetic term of 
our action becomes 



2 



1 ^TrQ^r; 1 - 0) 2 + im 2 Tr0 2 - 



\m* + £(1 - cos(^) 



<f> q 4>- q e- 2lti ^l N , (2.13) 



coinciding with the standard lattice propagator of a scalar field on a size L D- 
dimensional lattice up to the last phase factor in (2.13). 

The role of this and similar phase factors appearing in the vertices of the 
theory is to restore the topological expansion in large- N. We have now obtained 
an effective field theory with scalar variables <p q . The phase factors keep track of 
the fact that the theory came from a matrix field theory and that the expansion 
in powers of 1/7V can be represented in terms of a topological expansion where 
the leading order corresponds to spherical topologies, and higher orders come from 
surfaces of increasing Euler number. 

In particular, the phases of all planar graphs vanish. The easiest way to see that 
is to note that the phase factor does not change when we contract all the propaga- 
tors connecting two vertices. The two vertices contribute factors Tr[A(q{) . . . A(q n )A(q[) . . . A(q' n ,)]^ 
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and Tr[A(— q[) . . . A(—q' n ,)A(q'() . . . A(q^,,)} while the phase factor coming from the 
propagators is exp J27=i(K\K)^j- ^ n together, the phase factor is 



5 (E * + E E «i + E ex p 



2ni 



+EW+E<^> - 



E 



E(^*i^') + E^^i 



i<j 



oc exp 



2ni 



where we used the 5-function constraints. This is precisely the phase factor of a 
single vertex with {fo}, {k'j} propagators attached. Provided we are dealing with 
a planar graph, we can repeat this procedure until we end up with a single vertex 
with no propagators attached whose phase factor is trivially zero. 

In case of non-planar graphs we will not be able to contract the propagators 
winding around non-trivial homology cycles. For example, it is easy to see that 
the phase corresponding to the torus diagram consisting of two propagators with 
momenta q\ and q 2 , each following one of the two non-trivial homology cycles and 
meeting in one vertex, is 

~2ni 



exp 



N 



({h\k 2 ) - {k 2 \h)) 



The non- vanishing phases of non-planar diagrams result in their l/iV 2 -suppression. 
Finally, note that (2.13) shows that the space-time degrees of freedom of the orig- 
inal theory are coded in the Fourier transforms of the "internal" SU (N) indices. 

The construction we have just carried out works only for even dimensions. 
For odd-dimensions we can leave one of the dimensions unreduced, and apply the 
reduction prescription to the remaining even number of dimensions. This yields 
an action 



S = 



dt^Tr<j)(tf 



dtTrV((f)(t)). 



We can consider a slightly more general model by including a general coupling in 
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the hopping term: 



±m 2 4> 2 -TrV Q {4>), 



Vb(0) 




m 2 = 2DR 2 + ml. 



(2.14) 



The hopping term is the first term in S((f>). Adding the effect of world-sheet 
curvature is easy, we simply replace the hopping term according to 



Our proposal to describe strings in dimensions higher than one is to investigate 
the properties of (2.14) and (2.15). We have effectively reduced the problem to a 
one-matrix model. The difficulty lies however in the presence of the twist matrices. 
We begin our analysis in the following section. 



Some general features of (2.14) and (2.15) are easy to extract. If the hopping 
term is very small, we are in the limit where the space-time lattice points are very 
far apart and the whole surface is mapped into a single point, or various points 
incoherently. This is the pure gravity regime, 7 s t is expected to be —1/2 and the 
induced metric plays no role. As the hopping term is increased, it becomes more 
likely for the surface to occupy more and more space-time lattice sites, the induced 
metric and singular embeddings in principle begin to play important roles, and 
we can expect a transition to a phase different from pure gravity at some critical 
value of the hopping coupling R 2 . We will show presently that in the most naive 
approximation we are driven to a branched polymer phase with 7 s t = 1/2, and 
at the transition point 7 s t = 1/3. A more careful analysis will show a much 
broader spectrum of possibilities. Note that the role played by R is similar to the 
compactification radius in the study of c = 1 strings propagating in a circle of 
radius R [29]. 



Tr(I>rt - 0) 2 h- TrAiT^Tl - 0)A(I>r+ - 0) 



(2.15) 



3. PROPERTIES OF THE REDUCED STRING ACTION 
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To deal with (2.14) we decompose the matrix <fi into its eigenvalues and angular 
variables as usual 



<p = U 1 \U A = diag(Ai, A2, • • • , Ajv), 

The angular integration produces an effective action 

e r cff «0 = J dUe R 2 E,Tr\r,(u)\rHu) 

,r eB [j>]-v (<t>) 



The simplest representation of F e g is given in terms of the leading order expansion 
in R 2 : 

Z = Jd<Pe N ( R2 ^ T ^) 2 -V^)). (3.2) 

The critical properties of (3.2) are very similar to those of the curvature model in 
[17]. What the first term in T e g is describing is the breaking of the surface. The 
Feynman graph representation of a term like (Tr0 n ) 2 is a vertex where two surfaces 
osculate at one point, and the coordination number of the touching vertices in each 
of the surfaces is n. If we have instead a term of the form (Tr<f) n ) p , it represents 
p-surfaces sharing one point. When there are no products of traces, the Feynman 
graph expansion for the free energy or the connected correlators involves only the 
sum over connected orientable surfaces. When products of traces are included we 
obtain an expansion where we have surfaces of different topologies touching at 
points determined by the new vertices. The new graphs look like "daisies" . In the 
original curvature model introduced in [17] the initial partition function is 

Z = J d^-N^A+gTr^) ^ 

The angle integral is again very difficult to carry out in this model, and the simplest 
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non-trivial approximation is to study the action: 

S cS = N Qtt-0 2 + gTrtj) 4 + ^(Tr0 2 ) 2 ) . 

For real A, the coupling g' = (TrA 2 /N — (TrA) 2 /N 2 ) is always positive. In this 
case as we will see later one cannot escape the pure gravity phase. To leading order 
in 1/N, the graphs have the topology of spheres touching at one point, where any 
two spheres in the graph can share at most one common point, and no closed loop 
of spheres can appear. A generic graph looks like a tree of spheres. In the planar 
approximation it was shown in [17] that we have the following possibilities: 

1) . g' > 0, 7 s t = — 1/2 and we are in the pure gravity regime. 

2) . g c < g' < 0. Again 7 s t = —1/2 and we continue in the pure gravity case. 
The branching is still unimportant, and the leading contributions to the free energy 
come from graphs with few big spheres and few small branchings. 

3) . At g' = g c = —9/256, 7 s t = +1/3 and the number of branchings become 
competitive with the number of smooth spheres. 

4) . For g < g c we obtain 7 st = +1/2 and the dominating term is the one with 
a product of two traces, thus the notion of a smooth surface disappears and we 
end up with a branch polymer phase. 

The same analysis can be carried out with our naive approximation (3.2), and 
we will only note the differences with respect to the analysis in [17] we have just 
reviewed. The basic differences are: 

1) . Our potential is cubic, Vq = ^0 2 + g<f) 3 /\/N, the one-cut solution is asym- 
metrical, and the new vertex is just given by (Tr0) 2 . These are purely technical 
differences and nothing in the analysis changes substantially. 

2) . We are always in the g' < region. In the original model in [17] we can 
reach this region only after analytically continuing the matrix A i— > iA in order to 
obtain the polymer phase. However, in this case the interpretation of the matrix 
A in terms of local world-sheet curvature loses its meaning. 
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Our simplified model then contains a pure gravity phase, an intermediate phase 
with 7 s t = 1/3, and a branched polymer phase with 7 s t = 1/2. We should like to 
remark that the exact solution to the pure curvature model (3.3) is known only 
for Penner's potential [30], however in this case the full theory has the symmetry 
A i— > const. A -1 making the interpretation of A in terms of world-sheet curvature 
rather doubtful. 



4. PRELIMINARY ANALYSIS OF THE EFFECTIVE ACTION 



We have learned that the basic problem is the evaluation of the angular integral 



The naive approximation embodied in (3.2) already gave us a polymer phase and 
a critical point in-between pure gravity and polymers. To understand the behavior 
of strings in dimensions higher than one, we need to obtain information about the 
general critical properties of the effective action in (4.1). Since R 2 is a hopping 
parameter, we can write the hopping term as = (f)T ^(pT^ and in analogy with 
standard lattice analysis, we could write an approximation to the effective action in 
terms of a sum over connected graphs of angular averages over the hopping terms 
generating the graph: 



where the subscript U indicates that the average is taken only over the angular 
variables. To every connected graph the average over angles associates a set of 
products of traces. We can interpret these traces as different ways of fracturing 
and breaking the loop C. They describe how the embedded surface is collapsed 
into pieces. If we add a local curvature term (the A- matrix), we will obtain more 
general operators. Since we do not have an explicit way of evaluating (4.1) and 




(4.1) 



r eff = (TrY[L(l))u, 



C,Connected IsC 
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since we are interested in the critical behavior of (2.14) and not in the fine details of 
lattice dynamics, it is worth exploring the type of universality classes of potentials 
with arbitrary numbers of products of traces. There are two qualitative classes of 
potentials, depending on whether they contain a finite or infinite number of traces. 
So far we have analyzed the case of finite numbers of traces. 



Defining 



U = ^Tr<t> 2 \ (4.2) 



we take the effective action to be a general function of an arbitrary, but finite 
number of ti variables, 



r eff = V(h,t 2 , ...,t n ) = V^Tr^ 2 , lTr0 4 , . . . ^Tr<p 2n ). (4.3) 



Thus we are faced with the analysis of the critical properties in the planar limit of 
the action 

S(<j ) ) = N 2 V(t i ) + NTrV (<j ) ), 
k 

The world-sheet cosmological constant is related to the the coupling gi = e~^ B . 
We have made the simplifying assumption of restricting our considerations to the 
case of even potentials. This is just a question of technical simplicity, and nothing 
changes if we take general potentials. The analysis can be carried out for any of 
the Kazakov potentials [31]. 

The solution of the planar approximation of (4.4) uses the Hartree-Fock (HF) 
approximation, which becomes exact in the planar limit. We first write down 
the saddle-point equations for (4.4), and then we replace the traces (4.2) by their 
planar vacuum expectation values X{ = {jjTrcj) 2 ' 1 ). This reduces the problem to 
the original pure gravity case studied in [24]. Finally we fix the variables Xi self- 
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consistently. The saddle-point equations and the HF conditions are given by 



-^2^ p A 2fc - 1 = P.V. / dii 



p(p) 



X — ii 



9 P = 9 P + 



dV(x) 
dx k 



(4.5) 



= J P(/V fe - (4-6) 



p(/x) is the density of eigenvalues. We will denote by V the potential with the 
shifted couplings (4.5). 

Following the pure gravity case, we look for one-cut solutions of (4.5). This 
means that the loop operator 



F(p) = Jdfi 

is given by 



Pip) 
p-fi 



F(P) = \v\p) ~ M(pWp>-R. (4.7) 

The polynomial M{p) is determined by requiring that at large values of F{p) ~ 
1/p + 0{l/p^). Since V depends only on a finite number of traces, the potential 
V is a polynomial, and there are no ambiguities in the determination of M. 

Define 

G(t) = P F(p)\ p -2 =t = 1 + Xl (R)t + x 2 {r)t 2 + ■■■. 
The condition xo = 1 is the string equation. Since G(t) is analytic around t — 0, 

j>t l G{t){l-tR)- 1 / 2 = 0, (4.8) 
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for I > and C a circle around t — with radius less then 1/R. Using (4.7) and 

the equation (4.8) results in M(p) = z~2 m >o M m p m with 

Mm = Yl ( J ) ( m + ^ + 1) W+i- 

fc>0 ^ ' ^ ' 

The density of eigenvalues is given by 



(4.10) 



p(\) = 1-M(\)\Jr - A 2 , A e [->/!, >/!]. (4.11) 
In 

In the following we will also need expressions for x p (R) = (-^Tr0 2p ). To derive 
those, we find it convenient to distinguish between two sources of -R-dependence: 
the explicit independence that would be here even for pure gravity, i.e. for = g^, 
and the implicit independence of x p (R) (and therefore of gk), generated by the HF 
feedback. Consequently, we will use two kinds of derivatives w.r.t. R: the usual 
derivative d/dR = Or, taking into account both sources of -R-dependence, and Dr, 
derivative w.r.t. R carried out while ignoring the implicit HF dependence of x^. 
In particular, 



D R G(t) 



VI -Rt 



2 " v ' ' \ t 



(4.12) 



Since yl — RtDnG(t) is a series in positive powers of t, the only term from the 
square brackets on the r.h.s. of (4.12) that contributes is Mq/2 + RDrMq. Thus, 



Drx p = <j> 



Vl - Rt ' 

^™)( 2 ;)(i) p ' 
2 ;) (?)'»- 
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Finally, using the explicit expression for Mq from (4.10) we arrive at 



The string equation is given by (4.14) with p — 0. We are interested in the 
analysis of the critical points of (4.4), the computation of string susceptibilities, 
and in the computation of correlation functions of macroscopic loop operators. If 
we denote by (3 the coupling g2, the string susceptibility is determined from 

dx2 0x2 dRdx2 
X = ~df3 = W + d(3~dR' ' 

This is because the coupling (3 multiplies the operator £2 i n the action. Therefore 
dh\Z/d(3 is proportional to the expectation value X2, and the second derivative of 
the logarithm of the partition function is the string susceptibility. 

The dependence of R on (3 is read off from the string equation xq — 1. In 
particular, 

dxQ dxQ dRdxQ 
= ~dj = ~d(3 + 'df3~dR 

allows us to read the equation determining the critical points, 

d(3 d R x 



dR dpXQ ' 



(4.16) 



Since the denominator in (4.16) is proportional to R 2 , the critical points are de- 
termined by the zeroes of the numerator. We can similarly write \ according to 

X = dpx 2 - dpXQ |^ . (4.17) 

The singularities in the behavior of \ come from the ratio of derivatives in 
(4.17). For the standard Kazakov critical points [31] (4.17) does not blow up at 
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the critical point R c , and x ~ Xc ~ R — Rc- Since at the m-th critical point 
(3 — f3 c ~ (R — R c ) m this yields 7 s t = — 1/m. As we will show presently, the critical 
points we describe have 7 s t > and x will be singular at the critical point. All the 
phase transitions are second order. 

To derive the criticality conditions, use (4.14) to write the complete derivative 
of x n w.r.t. R as 




(4.18) 



Defining the matrix 




(4.19) 



we can compute d R x p in terms of D r xq: 

d R x p = [(l-U")- 1 A] p D R x , (4.20) 

where A is a vector whose p-th component equals (^) When p = 0, we 

obtain the criticality condition implied by (4.16): 

d R x = [(1 - U'^A^Drxo = 0. (4.21) 

Notice that the criticality condition splits into two terms. The first one is 
related to the function including the couplings with multiple traces, 

[(1 - U"Y l A\ = 0, (4.22) 

and the second is equivalent to the criticality condition for the Kazakov critical 
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points: 



Drx = 0. 



(4.23) 



The term in (4.17) which may lead to singularities in x an d to positive 7 s t is given 
by 

d R x 2 [(1-U")- 1 A} 2 



d R x \{l-U")^A] Q 



(4.24) 



In (4.21) there are three possibilities: 1). The polymer couplings in U become 
critical but the gravity part does not. The zeroes of d(3/dR come only from the 
polymer contribution. 2). The gravity contribution is the only one generating 
zeroes of df3/dR. 3). Both terms become critical. In the first case we have a theory 
of polymers, and little is remembered of the coupling to gravity. In the second case 
the polymer degrees of freedom are completely frozen and we reproduce the string 
susceptibility and exponents of the Kazakov critical points. The third and more 
interesting case is when both the "polymer" matter and gravity become critical 
simultaneously. This is the case where we can find novel behavior. Notice also 
that unless we tune the numerator in (4.24) we will generically obtain 7 s t > in 
cases 1) and 3). 

If R c , p c are the critical values of R and (3, and if near the critical point 

[(1 - U")- l A] Q ~ (R - R c ) n D R x Q ~ (R - R c ) m , (4.25) 

then 

P-P c ~{R-R c ) n + m + 1 . (4.26) 

If we do not tune the numerator of (4.24) to partially cancel the zeroes in the 
denominator, we obtain 

1 

X ~ (R-R c ) n ' 
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The previous two equation give a string susceptibility for the (n, m)-critical point: 



7st = —^—7- (4-27) 
n + m + 1 

By tuning the numerator in (4.24) we can change the numerator in (4.27) to any 
positive integer smaller than n 

P 

7st 



n + m + 1 

We will show below that only in the case n — 1 the model maintains some 
resemblance with the properties one would expect of a non-critical string. We 
can define operators creating macroscopic loop only in that case. For n > 1 
the polymerization of the surface is so strong that there is no room left to open 
macroscopic loops. 

To illustrate the preceeding discussion, we present an explicit example. The 
simplest case to study is the one where the function V depends on a single trace 
V = V(ti). In this case the matrix U" becomes: 

The matrix U" contains a single non- vanishing column in position I. The inverse 
of 1 - U" is: 

U" 



(1 _ U")~ = 1 + , 

i-WTM!) 1 ) 2 ^ 



therefore, 



\{l-U")-'A] p = 

2 

'2p\{R\ p 1 I 2 ( (2l\ [R\ l \ (2p\(R> l> 



V/ i-|(( 2 /) (%)yv»P + l \\ l J V4y / \i> 
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l+l(( 2 /) (f)>" 



[(l-U")- 1 A]o = 



i-K( 2 /) (f)W"' 



(4.28) 



The criticality condition (4.22) is now 

2 / r>\ 21 



(4.29) 



We may simplify the expressions for the moments by eliminating V' from the string 
equation xq — 1 to get 



l + k\k 



2k 



pk(p-l) (2p\ (R 
p + k \p 



9 P 



(4.30) 



The analogue of the criticality condition (4.23) is therefore 

rdrxq = i + J2p(p-i) ( 2 j) (f )%p « (fl - ^c) ? 



(4.31) 



Note that we still need the string equation for / 7^ 2, since (3 = gi appears on the 
right hand side of these expressions. Hence, the HF feedback is explicitly solvable 
only in the case I = 2. On the other hand, direct differentiation of (4.30) gives 

dxk k 



dR l + k 



A k D R x . 



(4.32) 



This means that the total derivative dxi~/dR is proportional to the critical- 
ity condition (4.31), thus restricting the critical points accessible with the simple 
potentials V(ti) to (n, m = 0) or (n = l,m). The criticality condition (4.29) for 
n = 2 implies 

2 / ,\ 2 



dW 



dR 



= V 



R — R c 



( 2 /)(^n^)>4(( 



Rr 



dx\ 



dR 



R — Rc 



-21 + 



I 



2 A ( R, 



V"\ Xl (R c ))D R x Q (R c ) = 0, 



where W = | ^( 2 ^) (-j) j V"(xi(R)). Obviously, the last equation is inconsistent 
with Drxo(R c ) = 0, the condition for m > 0. One can show, following the same 
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line of reasoning that we used here, that a more general class of potentials with 
V = V(^2 a r t r ) can accomodate n > 1, m > points, with a judicious choice of 
a r 's. Since the physically interesting case will anyhow turn out to be n — 1, we 
will relegate the discussion of this more general class of potentials to appendix A. 

Finally, let us display a few explicit potentials resulting in a critical behaviour 
of the new kind that we are studying. Specializing to V = Vfa), the criticality 
conditions (4.22) and (4.23) become 




(4.33) 



where R — R c (l + x) and x — ■> is the critical point. Another useful identity is 

-»(s ,+ »g^(?)(r- <-> 

It will be sufficient to look at the potential 

U = g 1 t 1 +pt 2 + htl 

We choose R c = 4. To obtain n — 1, m — critical point all we need now is 
h = -1/72. The string equation x = 1 fixes /? = 1/18 + x 2 /Q + 0(x 3 ) which 
implies 7 = 1/2. The eigenvalue density (4.11) is (neglecting 0(x 2 ) in M(A)), 

^)=K5-f-f A V 4(i+i) - AS ' 

which for x — > reduces precisely to Wigner's semicircle law . For n — m—1 (4.33) 
gives /i = -1/72 and g 1 = 1, whereas x = 1 gives /3 = -1/18 - 2x 3 /9 + C(x 4 ) 
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implying 7 = 1/3. The critical eigenvalue density is exactly equal to the pure 
gravity one (7 = —1/2) 

P(A) = ^(4-A 2 ) 3/2 - 

Note that the critical eigenvalue distribution reflects only the Kazakov part of the 
criticality conditions, while it is insensitive to (4.22), the condition responsible for 
divergence of x an d for positive string susceptibility. 

An argument similar to the one presented at the end of [17] shows that the 
critical point for the transition between the polymer and the pure gravity phases 
is generically second order. To see this, we look at the partition function for 
"daisies" at a fixed number of Feynmann vertices (or elementary plaquettes in the 
dual version). We can extract this quantity from the asymptotics 




So that, in the thermodinamic limit n — > 00 we have the free energy 

/ = log|&| 

where (3 c (g, V) is the convergence radius of the planar expansion as a function of 
the couplings in the matrix-model Lagrangian. g denotes collectively the linear 
couplings except #2 = P, while V corresponds to non-linear terms like (Tr0 2 ) 2 . 
We are interested in singularities of /3 C as a function of V, since these couplings 
act as the "temperature" for the gas of splitting points. A first order phase tran- 
sition would imply coexistence of the smooth and polymer phases, while a second 
order transition represents global condensation of daisies at the critical point. We 
calculate the function f3 c (g, V) as 

P c (g,V)=P(R c (g,V)),g,V) (4.35) 

where (3(R,g, V) is the solution to the s self-consistency conditions (equations for 
the moments Xk) plus the string equation. It depends on the planar solution 
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and it is a regular function of the couplings g, V. R c (g, V) is the position of the 
critical point as a function of the couplings, and we can find it from the criticality 
conditions (4.21): 

^(R c , g ,V) = (4.36) 

Tunning g and V to some submanifold of appropriate codimension we can arrange 
zeroes or any order in (4.36). Further change of g, V within these universal sub- 
manifolds amounts simply to move the location of the critical point R c - In our 
case (4.36) has a factorized structure as in (4.21), and we have two branches of 
critical points: R\ and R]} , determining two functions in (4.35); fi\ and (3^ . The 
smaller of these values determines the radius of the covergence disk for the planar 
perturbation series around (5 — 0. Hence we are interested in the behaviour of the 
function 

[3 c (g,V)=mm{\[3 I c l\p I c I \} 

near the coalescing region Rj, ~ R]J ~ R*, reached by tunning V. 

(3 c (g, V) is continous across the transition, and the first derivative is 

§ to. v) = % (W) + f§W 9 , vU,v)^( a , v) 

The first term on the right hand side is continous. In the second term, may 
have some insolated singularities, but we can assume smooth behaviour in open 
sets. On the other hand f^(-R c ) vanishes identically for all critical points, not only 
for the coalescing one. So the second term is zero along the critical submanifold 
and is continous at V* . 

The second derivative contains unprotected terms like q^q V (Rc(g, V), g, V). 
These are generically discontinous since R^ depends on the first derivative of the 
potential, while R]J depends on the second derivative. Thus we are led to at least 
second order phase transitions and a continuum limit is possible as expected. 
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We have presented the general analysis of effective actions in the planar limit 
containing arbitrary, but finite numbers of traces (4.4). We were able to show that 
beyond some critical coupling, 7 st becomes positive and takes values of the form 
7st = n/(n + m + 1) for any positive integers n, m. We should continue exploring 
the properties of these critical points and their scaling operators. Obviously the 
more challenging case of studying arbitrary functions V depending on an infinite 
number of traces is crucial before we can draw any conclusions concerning the 
properties of our model. This together with the study of general properties of (4.1) 
should shed further light into the properties of non-critical strings in interesting 
dimensions D = 2,3,4. It is very important to obtain in what way the effective 
actions studied in the previous sections depend on dimensionality. 

5. MACROSCOPIC LOOPS 

To explore the new critical points in some detail, we compute the correlation 
functions of macroscopic loop operators. The simplest loops to compute are the 
ones analogous to those appearing in pure gravity (see for instance [32]). They are 
obtained by taking the limit k — > oo of Tr<fi 2k . In terms of the original reduced 
model these are loop on the surface whose position is averaged over the target 
space. One could also define loops with definite positions in the target space. The 
expectation values of the latter would certainly give crucial information on the 
reduced formulation of string theory, but we have not yet calculated them. For the 
simpler, pure gravity loops of length / we take the definition 

(w(l)} = lim ,/^(Tr0 2fc ) = N lim J^x k . (5.1) 

As a matter of convenience we choose the critical values for f3 and R as f3 c = 
— 1,R C = 1. In the Kazakov critical points, 7 s t = — 1/m and the scaling limit is 
taken to be iVa 2_7st = k~ 1 . It is a simple but non-trivial fact that the correct 
way to define a loop of length / in this case is to scale k according to ka? lst = I as 
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k — > oo (see for example [32] and references therein). In our case we will restrict 
for simplicity our computations to the case of potentials V = V(ti). 

We will start by calculating dx^/dgj. From (4.14) we have 
dxi _ j+i( 2 /) Cj) (f Y 



d 9j l-W 

d_x, =J ^ (2k\ (2 A (RV +k 

dgj j + k\k)\j)\Aj 1-W 



(5.2) 



where we have introduced the notation W = ^ ^( 2 ^) (^) j V"(xj). From (4.28) 
we have 

d R x k _ (2k\ (R\ k 

Putting it all together, 

cbfc _ dxk_ _ drrp (9j?x fe 
dgj dgj dgj d R x ' 

jk (2k\ (2 A /iA jk (2k\ (2 A /iA 2/H/ 

(5.4) 

For V"" = 0, up to irrelevant numerical factors this answer coincides with the 
one-loop expectation value in the Kazakov critical points. 

To go to the continuum limit we have to choose the scaling variables. Recall 
that in the double scaling limit [11] the free energy takes the form 

where \x = — log j3, and g is the genus of the triangulations summed over. The 
renormalized cosmological constant is introduced according to /i — /jl c — a 2 t where 
a is the microscopic lattice spacing. Then the effective string coupling constant is 
defined by Na 2 ~ lst = k~ x . Hence the scaling variable associated with the string 



28 



susceptibility is x = a 27st «- At the (n, m)-critical point x behaves as x ~ (1 — 
R)~ n . Thus the scaling variables are 

X ~ t ^r- = a~ 2 ^u, 1-R = a 2 W^-iM 

(1 — it!) 71 

/3 - /3 C = /3 + 1 = a 2 t. (5.5) 
For an (n, m)-critical point we have 

2 — 1 In r> n -i 

1 — i? = a n + m + l u n , x = ua n+m+l , Na »+m+i = ^ . (5.6) 

An insertion of a loop operator is a string interaction which brings in an extra 
power of k, the string coupling constant. Therefore, in general we expect 

(w{h) . . . w(£ M ))c = K M - 2 {w(ti) . . . w(£ M ))c, (5.7) 

and in particular 

(w(£)) = -(w(£)). (5.8) 

K 

This will allow us to arrive at the proper definition of the length stick, i.e. to fix y 
in£ = kay. Setting j = 2 in (5.4) (g 2 = (3, 1 + W oc (R c - R) n ), we easily obtain, 
up to a non-universal positive multiplicative constant 

4(Tr0 2fc ) = Na 2 (^l(l- a^nu^Y^ uaT^ | . (5.9) 
Identifying the r.h.s. with \^[{w{£)) leads to 

2n 

£ = ka "+m+i ; 

4>(*)> = -^-e- £ / M forn = l, (5.10) 
= for n > 1. 

In the phases with n > 1 the polymer couplings dominate completely the critical 
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behavior, in spite of the fact that gravity becomes also critical, in such a way that 
there is not enough room to open macroscopic loops* 

In the case of n = 1 we can also calculate two-loop and multiloop correlation 

2n 2n 2 _, 

functions. Let i\ = ka n + m + 1 , £2 = ja n + m + 1 . For n = 1, 1 + W = 2la m + 2 u . For 
k,j^oo and a — > 0, (5.4) gives the continuum two-loop operator 



M*iM* 2 )>= lim -^=(Tr^ k Tr^ 



e -(^i+^)/« 



?1 + ^2 ^2 



^_ e -(^l+^2)/«_ 



(5.11) 



The lesson we can draw from (5.11) is that in the simple case of n — 1 the effect 
of the polymer couplings is to contribute an extra state (the last term in (5.11) ) 
which resembles very much the contribution one would expect of a tachyon. This 
term represents the breaking of the cylinder interpolating between the two loop of 
lengths li, I2 into two disks osculating at one point. We can represent the last term 
in (5.11) as 

Mh)P)jj^(Pw(l 2 )) 

where P is the puncture operator. 

With (5.10) and (5.11) it is easy to obtain the multi-loop correlators as well. 
We can write {w(£i)w(£2}) in two ways: 

00 

MW»» - fa >- (-£) ^ / (l + f) ^-"<"< (5-12) 



and 



{w{h)w{t 2 )) 



4/K)R)§- ,f ' +fe,/ * 



* Strictly speaking, we should have put this discussion in the framework of potentials with 
V(5^a r t r ), the simplest ones with n > 1 critical points. As shown in the appendix A, the 
much more complicated formulas in that case lead again to (5.9) and to the same conclusion 
that (w(£)) = for n > 1. 
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Comparing (5.12) and (5.13) we have 



lim 




Define the operator of insertion of a macroscopic loop D{£) as 



D 




d 



Where = As it should be, D(£) is proportional to k, and we can write 
the M-loop amplitude as 



The pattern is reminiscent of the residue of the tachyon pole in the Belavin- 
Knizhnik theorem [33]. This interpretation, though tempting, has to be taken 
with several grains of salt because in the types of effective actions we are studying 
it is hard to see any dependence on dimensionality. Nevertheless we take this re- 
sult as encouraging evidence that our approximation captures some of the expected 
properties of non-critical strings for D > 1. 

The form of the two-loop operator (5.11) does not depend on the form of the 
potential V(ti, . . . , t n ) as long as we consider (1, m)-critical points. Hence only the 
quadratic part of V determines the critical properties of macroscopic loops. 

To conclude this section we mention that we can study, at least in part, the 
spectrum of the scaling operators by perturbing the criticality conditions. The 
string equation is then modified according to 



The subleading contributions depend on both the polymer and Kazakov perturba- 
tions of the critical conditions (4.21). The bare couplings of the scaling operators 
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! + /? = (!- R) n+m+1 + T ki l ~ R ) 



k+l 



are the 77? parameters. Using (5.6) the string equation becomes 

n+m + l V -V fc+1 . . 

t = u « + y. T k u n , (5-14) 
k 

and the 77 are the renormalized couplings. In the original model these couplings 
would hardly exhaust the scaling operators, although they may represent a signif- 
icant subset associated to the simplest loop operators. If we call the scaling 
operator associated to the coupling 77, its planar correlators at the (n, m)-critical 
point follow from (5.14) 

du _ fc+i du 

= — U n — 

dTk dt 



and 



Tb k — m — 2n 

(a^PP) = 1 n+m+l . 

n + m + 1 



6. l/7V 2 -CORRECTIONS, DISCRETE 
AND CONTINUUM LOOP EQUATIONS 

In this section we will use Schwinger-Dyson equations to evaluate l/iV 2 -corrections| 
to our planar results. The l/iV 2 -corrections will prove to be consistent with string 
perturbation expansion and with double scaling* In particular, we will calculate 
the string susceptibility on the torus to be 

7torus = +2. 

It is amusing to note that this is in agreement with DDK formula 7^ = 2 + (1 — h)jQ 
for string susceptibility at genus h, where 70 (a complex number for c > 1) drops 
out at h — 1. 

* Since the equivalence of the original action (2.3) and its reduced version (2.12) has been 
demonstrated only in the planar approximation, the l/iV 2 -corrections for (4.4) are one step 
further removed from (2.3). Still, in order to be able to take the new critical points of our 
effective action (4.4) more seriously, we want to see if their continuum limit is consistent 
with a string theory. 
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We will study matrix integrals of the form 



/ 



Z = / £>0exp ^ -N z 



.n=l ^ 7 



1 ( 1 V 

m=2 ^ 7 



The structure of the calculation follows [34]. In a standard fashion, 



0= / P0^j(e^Uexp<!-iV 2 



leads to the loop equation 



.n=l ^ ' 



* / / 1 \ m_1 1 

m=2 \ v 7 

Let us define connected correlation functions ^£ 1 ...g M as 



Tr 6 6 e l 



(6.1) 



= N M - 2 ( Tre 1 ^ ■ ■ -Tre tM ' 



where the explicit factor of iV makes all <^£ 1 „4 M of order one. For N large, (6.1) 
leads to 



Tl — 1 n 



+ 



AT2 



- 2V"'{x 2 )^ A l^Tr (*V*)\1 + (9(l/iV 4 ), 



(6.2) 
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where 



9n = 9n + V'(x2)5 n ,2, 
1 



N 



Tr<p 



2i 



$ n , m = {Tr<j ) n Tr<f ) m ) l 



After the Laplace transform, (6.2) becomes 



V'(p)G(p) = Q(p) + G 2 (p) + 



N 2 



$(p,p)-4y"(x 2 )( P 3 $ 4 (p)-$2,4) 



2V'"(x 2 )^Ap 3 G(p) - V 2 - xi) 



+ 0(1/N 4 ). 



We have introduced the notations 



n=l 
s n—\ 



Q(p) = ^n~g n J2p 2n ~ 2t ~ 2 ^ 



oo 



Tr- 





oo oo 







$ 4 (p) = (Tr^Tr^—) = $4,2nP _2n_1 - 



The solution to (6.3) is well-known in the planar limit: 



G {0) ( P ) = \v\p) - M(p)y/p^~R, 



n=l n=l fc=0 V K J V J 
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where R is fixed by the string equation 



x = 1 = ^2ng n 



'2n\ fR xV 



n ) V 4 



To include l/iV 2 -corrections, write 

V> = V'(°) + SV' = V'M + AV"(x 2 )5x 2 p 3 , 
G = G (0) + SG, 
Q = Q(°) + 5Q, 

where 5( ) is 0(1/N 2 ). From (6.3) we obtain the equation for 1 /^-corrections 
4V"(x 2 )5x 2 p 3 G^(p) + 2M(°)(p)v / P 2 - R5G(p) = SQ{p)+ 



1 

+ iV2 



(0) 

$(p,p) - 4\/"(x 2 )(p 3 $ 4 (p) - $2,4) - 2V w "(:r 2 )$4,4(p 3 G(p) - p 2 - x x ) 



(6.4) 



where 

s n— 1 

<*Q(p) = E 2n &» EP 2n_2i_2fai + 4^"(* 2 )MP 2 + xi). 



n=2 



An important piece of information that we can extract from (6.4) is the de- 
pendence of 5x 2 on R — R c . This will give l/iV 2 -corrections to specific heat. To 
simplify (6.4) we use the fact that, for s > 2, M^°\p) vanishes for p equal to some 
Pq. Therefore, at p = po (6.4) reduces to 



4V"(x 2 )5x 2 p 3 G^(p ) = 5Q(p ) + ± 



$(p ,Po) - W'(x 2 )(pl§±(po) - $2,4) 

(0) 



-2V"\x 2 )^^{plG{p Q )-pl-x 1 ) 

(6.5) 

We need the planar expressions for and $4(p). They are calculated in 
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appendix B. The results are 



$ 4 (p) 



1 + W 
R 



8p 4 - ARp 2 - R 2 n o 

; %P 

A^fp^R 



A(p 2 _ i?)2 
9i? 4 (l + W) 



6 16p 3 (i? 2 + 4V-8p 4 ) i? 3 (8p 2 + it!) 
izop H ; = h 



Vp^R 



p 2 -R 



(6.6) 



We are ready to look at some explicit examples. As discussed earlier, V(x 2 ) 
potentials can describe critical points with either m = or n = 1. The simplest 
non-trivial case are potentials with s = 2 which allow us to look at m = 0, n > 1 
and m — n — 1 points. 

For s = 2 (6.5) gives 

(PoG(po) - Po - zi) (V'(z 2 )^ 2 + l_V"\x 2 )^ = 
= A~g 2 5xi + ^ [$(po,Po) - 4\/ // (x2)(po$4(po) - $2,4)] • 

Another relation between the unknowns 8x\ and 5x 2 comes from the coefficient of 
1/p 2 in (6.4): 



Sx! = —(92 + V"(x 2 )x 2 )5x 2 - ±-^(2V"{x 2 ) + V'"(x 2 )x 2 ) 
91 N z gi 



Finally, 



A 1 M 
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where 



N = 



R 



9R 4 



A(p 2 - R) 2 16(1 + W) 



ZRfcV'M I G(po) h 

9M-R) + {2) { 2 + ~ 91 X2 



V = AV"(x 2 )G( Po ) + ^(~g 2 + V"(x 2 )x 2 ), 

91 

G(Po) = Po G (Po) ~ Po ~ x i = 9iPq + ^92Pq -Pq~x\, 



2 R 



91 
^92 



Specializing further to m = 0, n > 1 potentials with 70 = n/(n + 1), we set 
it! = 4(1 + x < 1 and #1 = 1/2. From xo = 1 we get #2 = -x/(12(l + x) 2 ) and 
Pq = 3/rr + 0(1). Note that #2 is 0(x) and Pq diverges as we approach the critical 
point, which is as it should be since m = potentials have Gaussian linear part of 
the potential (i.e. g n — 0, n > 1). 

Setting finally h,2 = —1/72, h m = 0, m > 2 we have fully specified the m = 0, 
n = 1 potential. We have M = -2x/3 + 0(x 2 ) and V = -Ax 3 /9 + 0(x 4 ). 
Furthermore, (3=1/18 + x 2 /Q + 0(x 3 ). Therefore, 



5x 2 oc 



1 1 



dSx 2 11 1 1 

OC ~-^—r OC 



d/3 



N 2 x 4 



N 2 ((3 - (3 C ) 2 



and 



71 = +2. 



We have derived the first two terms in the genus expansion of the specific heat. 
In terms of renormalized continuum variables we have 



d 2 



K 



h ft ~n + 

r?° t 2 



(6.7) 



The l/iV 2 -correction is consistent with string perturbation expansion and with the 
double scaling limit. 
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One can repeat this calculation for m = 0, n > 1 and one is always led to 
(6.7). Similarly, for n — m — 1 potential that we already discussed we have g\ — 1, 
h 2 = -1/72, /i m = 0, m > 2, g 2 = -i^fgjr, /5 = -1/18 - 2x 3 /9 + 0(x 4 ), 
dga oc ^ and we again have 7l = 2 and (6.7). 

To reach n — 1, m — 2 we need potentials with s = 3. Here M(p) will have 
two zeroes, (p^) 2 , that will both coalesce with R at x = 0. If we take g\ = 3/2, 
^3 = 1/60, h 2 = -1/72, h m = 0,m>2 we obtain /3 = -7/30 + x 4 /4 + C(a; 5 ) and 

(p+) 2 = 4 + (-1 + v^rr - 2^^ + 0(x 3 ), 
(p-) 2 = 4 + (-1 - v^rr + 2 V /Z 5a; 2 + 0(x 3 ). 

Now (6.4) will result in two equations with two unknowns, 5x\ and 5x 2 - Their 
solution again leads to 71 = 2 and (6.7). 

It is interesting that there is a string perturbation expansion (6.7) even for n > 
1. At those critical points, as we have seen through vanishing of the macroscopic 
loop operators, polymerization is so strong that it is impossible to open big holes 
in the surface. Since the only continuum objects used in the derivation of (6.7) arc 
microscopic ones like d/dt ~ P, we were able to avoid pathologies due to vanishing 
continuum macroscopic loop operators. 

Finally, let us look at the continuum limit of the loop equation (6.3) itself. 
Comparing the solution of the continuum limit of the loop equations with the 
continuum limit of the solution of the discrete loop equations is a useful check of 
our calculations. 

To begin with, use £ = ka n + m + l , p = p c + a n + m + 1 z, to determine the scaling for 
the universal, continuum loop correlation functions. On the one hand, we have 

(Tre^ . . . Tre lM ^j = k m ~ 2 I dz x ... dz M (w{ Zl ) . . . w{z M )) c . 
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On the other hand, 

T r Ji<t> Tre eM</) \ = - §e c = 

J- it .-.-Lie n M-2^Zi-£m 

= N M~2 J d V\--- dpM{w(pi) ■ ■ ■ w(p M ))c- 

Consistency demands 

2(2 m +n+2)-M(3n + 2m+2) 

. . .w(pm))c = a m+n+1 (iu(zi)...ii;(zm))c, 

with n > 1, m > and M > 0. In particular, 

2m — re + 2 

(universal part = a m+n+1 (w(z)) . (6.8) 

The negative sign of n in the exponent is the first sign of problems for n > 1. 
We start again with n — m — 1 point. Set 

/3 = /3 C + a \ 
R = R c {l + a 2 /H 1 / 3 ), 
g 2 = (3 + 2h 2 x 2 , 
p = p c + a 2 / s z, 
gi — 1 + ca t. 

Demanding that universal (non-analytic in t) pieces in x\ and x 2 scale as (P) fixes 
R c = 4, f3 c + 2/^2 = —1/12. The criticality condition for n = 1 gives /12 = —1/72, 
/3 C = —1/18. Finally, since we want to blow up the critical region between the 
branch points, we fix p 2 = R c — 4. 

Plugging all of this into solution of discrete planar loop equation we obtain 

<^)>|planar = + 2z)y/z-tW. (6.9) 

On the other hand, scaling the loop equation itself, we arrive at 

9((w(z)) 2 + k 2 (w(z)w(z)) c ) = (24 + 6c)t - 12(P)z + 16z 3 . 
With c = —14/3, (P)|pianar = t 2 ^ 3 we recover precisely (6.9). We have shown that 

39 



operations of taking the continuum limit and solving the loop equations commute. 

Now we turn to n = 2, m = 0. According to (6.8) G(p)| U niv. = (w(z)) and 
universal loop equation should appear at Oil) in scaled discrete equation. Indeed, 
for p 2 c = R c = 4 we obtain 



For n > 2, m = negative powers of a in (6.8) give immediately (w(z)) = 0. 



In this appendix we examine with some detail potentials of the form V(Y1 a r t r ) =| 
Ylim hm{oi ■ t) m . These potentials already contain complicated couplings between 
different traces, but they are still explicitly soluble and serve as multicritical po- 
tentials for m, n > 1 critical points. 

The general matrix in (4.19) simplifies in this case to 



= (w(z)} 2 + k 2 (w(z)w(z)} 

= (U>(^)) I planar- 



(6.10) 



7. Appendix A 



U' 



U <g> a 



where a, U denote the vectors of components a>k and 




(recall Ak = (^f) ) P ). This structure makes the inversion of 1 — U" trivial. Given 
that 



d^x = Drx + (U ■ a)djix = DrxqA + U(a • d^x) 



(7.1) 



we may solve d^Xk as 



dRX k = D R x Q 



(1-U- a)A k + U k a-A 
1-U-a 
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In particular 



9rXo = Dr ^° (i _ jj . a _|_ [/ oCt . ^) 
1 — U ■ a 



and the polymer criticality condition reads: 

p + A; 



1 - [/ • a + U a ■ A = 1 + V" J] -^-A p A fc a p a p ~ (R - R c ) n (7.2) 



Analogous considerations to those in (7.1) lead to 

24 fR\ 2 a-Dpx 



24 / R\ 
d P Xk =k^2 Ak \j) +Uk 



l-U -a 



With these ingredients at hand we can compute whose scaling for large k 
determines the one-loop correlator, 

dxk A Dgxo(l — U ■ a) + Uoa ■ Dgx 
d(3 k 1 + W 

for 1 + W = 1 — U-a + Uoa-A. Assuming an n-th order critical point for polymers 
in (7.2) we find the large k asymptotics 

dxk f%k\ ( R\ ^ji — RY n 



dp \ k J V 4 

the same scaling as for V(xi) potentials, and we conclude = here as well, 
for n > 1. 

Now we come to the definition of n > 1, m > 1 critical points. First we should 
simplify the expressions for the moments x k by solving V' from the string equation: 



V{a ■ x) = 1 Z^fVr 



and we get the functions 



x k = X k (R, g p , a p ) = ^ -^—^A p A k g p + 1 s^^J Jl^ A ^ a P ( 7 - 3 ) 

Notice that these functions do depend on g2 = (3, so that we still need the string 
equation to relate (5 and R. 
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The Kazakov-like critical points come from 

In this expression we assume gi = (3{R) solved from the string equation. 
Upon direct differentiation of (7.3) we get 

dX k n E P p%A p A k a p 

— — = D R x 

oR E P A p a p 

still proportional to Drxo, as in (4.32), but through a non-trivial factor. If R c is 
the critical point in the polymer phase then obviously 

dX k _dX k dX k d(3 dx k 

dfl " ~dR {Rc) + WdR {Rc) - ~dR {Rc) (7 - 5) 

Assume now a particular choice of a p so that, in the vicinity of R c = 1 we have 

Si? 22P A p( R c)a P 
for < r < m. Differentiating now (7.2) using (7.5) and (7.6) we arrive at: 
dW 



(R c ) ~ F(i? c )(l - i? c ) 2r + G(i? c )y'"(a • x)(l - R c ) 



dR 



with F and G regular functions at R c — 1. It is clear that dW/dR vanishes at 
i? c = 1 for < r < m. If r = m we need additional tunning of h m parameters in 
the potential so that V'"(R C = 1) = 0. Hence, we have explicitly shown how n = 2 
polymer critical points can coalesce with Kazakov critical points of any order m. 
In general, the possibility of having ^ L ^ at the critical point rules out the 
problem for arbitrary n. 

In order to find an n, m critical potential we are led to the following hierarchy 
of tunnings: we first tune a p as in (7.6). Then, for fixed a p , adjust g p to get (7.4) 
and finally we may tune h m in (7.2). It is interesting to note that the relative 
tunning among different traces in the non-linear potential was essential to achieve 
the full spectrum of string susceptibilities (4.24) quoted in the general discussion. 
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8. Appendix B 



In this appendix we will calculate and $4(p) needed in the calculation 

of l/iV 2 -corrections in chapter 6. The identity 



= / V<p^- i (e^) a6 lTre^exp <J -N< 



results in 

2n-l 



.71=1 ^ ' 



71=1 N 

+ E ^ (s) 3 O^ 4 ) = (8.1) 

J \N N N I N \N 



After expanding (8.1) in terms of connected correlation functions, using (6.2) and 
going to the planar limit one is left with 



V 7 <! (8.2) 



71=1 



The Laplace transform of (8.2) is 

2M(p)V^R$(p,p 2 ) = Q(p,p 2 ) + 9 G(p)-G( P2 ) 

OP2 P~P2 (8.3) 

- W"(x 2 )<Z> 4 (p 2 )(p 3 G(p) -p 2 - Xi), 

where 

s 2n-2 

Q(P,P2) = Y, 2n ~9n E P 2 ^ 2_ ^^(P2). 
n=2 i=l 
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For our even matrix integral measure, splits into even and odd parts. 



m,n—l 

oo 

= P~ 2m ~ 1 (pT 2n ~ 1 (Tr<f) 2m Tr ( i) 2n } c + 



m,n=l 



+ £ P" to (^)" to <rr^ 1 Tr^ 2n - 1 > c = $ev m (p,p') + $odd(p,p'). 



m,n=l 



We will calculate $ e ven using (5.4) to obtain (Tr4> 2 ^Tr(j) 2k ) c = —dx^/dgj for Z = 
2. To get $ oc id one can, following [34], expand (8.3). The coefficient of l/p 2 , 771 
will determine $2m-i(p)- Since all explicit dependence on V" drops out of this 
expansion, $ Q dd will have the same form as in [34]: 



_ , , 2p 2 R-R 2 

®odd{P,P) 



8p 2 (p 2 -R) 2 ' 
As for the even part, in calculating 



oo 



k,j=i d9 -i 

oo , 

dx k 



. y , 2*-l 

fc=l ^ 2 



we will need the following sums: 



8l ^ = 22kT2\k)\AJ ' 

fe=i v J 

»M=Tj+k{ j ){ k )ta) 
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One can evaluate them as follows: 




Now it is easy to show that 



$ 4 (p) 



1 



1 + W 
R 



lp 4 - 4Rp 2 -R 2 3 
4yV - R 



A(p 2 - R) 2 

AW 
9R 4 (1 + W) 



i oo 6 16p 3 (i? 2 + 4ify 2 - 8p 4 ) i? 3 (8p 2 + i?) 
128p H , = h 



p 2 -R 
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